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1. Introduction 

One of the important tools for studying the properties of quark-gluon plasma (QGP) in ultrarel- 
ativistic heavy ion collisions is the analysis of a QCD jet production. The medium-induced energy 
loss of energetic partons, "jet quenching", should be very different in the cold nuclear matter and 
QGP, resulting in many observable phenomena [Q]. Recent RHIC data on high-pj- particle pro- 
duction at y/s = 200 A GeV are in agreement with the jet quenching hypothesis (see, e.g., [§] and 
references therein). At LHC, a new regime of heavy ion physics will be reached at ^Ann = 5.5A 
TeV where hard and semi-hard particle production can stand out against the underlying soft events. 
The initial gluon densities in PbPb reactions at LHC are expected to be much higher than those at 
RHIC, implying a stronger partonic energy loss, observable in new channels. 

In the most of available Monte-Carlo heavy ion event generators the medium-induced partonic 
rescattering and energy loss are either ignored or implemented insufficiently. Thus, in order to 
analyze RHIC data on high-pr hadron production and test the sensitivity of LHC observables to 
the QGP formation, the development of adequate and fast Monte-Carlo tool to simulate the jet 
quenching is necessary. 

2. Physics model and simulation procedure 

The detailed description of physics model can be found in our recent paper [gj. The approach 
bases on an accumulating energy loss, the gluon radiation being associated with each parton scat- 
tering in the expanding medium and includes the interference effect using the modified radiation 
spectrum dE /dl as a function of decreasing temperature T. The basic kinetic integral equation for 
the energy loss AE as a function of initial energy E and path length L has the form 



where / is the current transverse coordinate of a parton, dP /dl is the scattering probability density, 
dE / dl is the energy loss per unit length, X is in-medium mean free path. The collisional loss in 
high-momentum transfer limit and radiative loss in BDMS approximation [Q] (with "dead-cone" 
generalization of the radiation spectrum for heavy quarks [§]) are using. We consider realistic 
nuclear geometry and treat the medium in nuclear overlapping zone as a boost-invariant longitu- 
dinally expanding quark-gluon fluid. The model parameters are the initial conditions for the QGP 
formation for central AuAu (PbPb) collisions at RHIC (LHC): the proper formation time To and 
the temperature Tq. For non-central collisions we suggest the proportionality of the initial energy 
density £o to the ratio of nuclear overlap function and transverse area of nuclear overlapping. The 
simple Gaussian parameterization of gluon angular distribution over the emission angle 6 with the 
typical angle of the coherent radiation 6q ~ 5° [Q| is used. 

The model was constructed as the Monte-Carlo event generator PYQUEN (PYthia QUENched) 
and is available via Internet [^J. The routine is implemented as a modification of the standard 
PYTHIA_6.4 jet event The event-by-event simulation procedure includes the generation of the 
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initial partem spectra with PYTHIA and production vertexes at given impact parameter, rescattering- 
by-rescattering simulation of the partem path length in a dense zone, radiative and collisional en- 
ergy loss per rescattering, final hadronization with the Lund string model for hard partons and 
in-medium emitted gluons. 

The full heavy ion event is simulated as a superposition of soft hydro-type state and hard 
multi-jets. The simple approximation [||] of hadronic liquid at "freeze-out" stage has been used to 
treat soft part of the event. Then the hard part of the event includes PYQUEN multi-jets generated 
according to the binomial distribution. The mean number of jets produced in AA events at a given 
impact parameter is a product of the number of binary NN sub-collisions and the integral cross 
section of hard process in pp collisions with the minimal transverse momentum transfer pf m . The 
extended in such a way model has been also constructed as the fast Monte-Carlo event generator [0]. 
Note that ideologically similar approximation has been developed in [JTot] . 



3. Jet quenching at RHIC 

In order to demonstrate the efficiency of the model, the jet quenching pattern in AuAu colli- 
sions at RHIC was considered. The PHOBOS data on T7 -spectra of charged hadrons [jll|] have been 
analyzed at first to fix the particle density in the mid-rapidity and the maximum longitudinal flow 
rapidity, Y™™ = 3.5 (figure [l]). The rest of the model parameters were obtained by fitting PHENIX 
data on p^-spectra of neutral pions [12] (figure ^|): the kinetic freeze-out temperature 7} = 100 



MeV, maximum transverse flow rapidity y™ ax = 1.25 and minimum transverse momentum of hard 
parton-parton scattering p™ m = 2.8 GeV/c. The nuclear modification of the hardest domain of pj- 
spectrum was used to extract initial QGP conditions: 7b = 500 MeV and To = 0.4 fm/c. Figure 
IH shows that our model well reproduces p T - and centrality dependences of nuclear modification 
factor Raa, which is defined as: 



Rhh{pT,r\;b) 



Og d^AA/dpTOlT] 

{Neon) d 2 a w /dp T dri 



where (N co \\) = 7\a(^) x a™ 1 is the average number of binary nucleon-nucleon collisions in a 
given impact parameter b (with nucleus overlap function 7aa(&))- If there are no nuclear effects, 
the value of 7?aa at high pj should be unity. 

Another important tool to verify jet quenching is two-particle azimuthal correlation function 
C(A<p) - the distribution over an azimuthal angle of high-/?/- hadrons in the event with 2 GeV/c < 
pj < p X j g relative to that for the hardest "trigger" particle with p T s > 4 GeV/c. Figure |] presents 
C(A<p) in pp and in central AuAu collisions (data from STAR Jl3|]). Clear peaks in pp collisions at 
A(p = and Acp = n indicate a typical dijet event topology. However, for central AuAu collisions 
the peak near % disappears. It can be interpreted as the observation of monojet events due to the 
absorption of one of the jets in a dense medium. Figure 4 demonstrates that measured suppression 
of azimuthal back-to-back correlations is well reproduced by our model. 

We leave beyond the scope of this paper the analysis of such important RHIC observables 
as the azimuthal anisotropy and particle ratios at low pj. In order to study them, a more careful 
treatment of soft particle production than our simple approach is needed (the detailed description 
of space-time structure of freeze-out region, resonance decays, etc.). 



3 



Simulation of jet quenching and high-pj particle production at RHIC and LHC 



LP Lokhtin 




Figure 1: The pseudorapidity distribution of 
charged hadrons in AuAu collisions for three cen- 
trality sets. The points are PHOBOS data, his- 
tograms are the model calculations. 



Figure 2: The transverse momentum distribution 
of neutral pions in AuAu collisions for three cen- 
trality sets. The points are PHENIX data, his- 
tograms are the model calculations. 
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Figure 3: The nuclear modification factor Raa for 
neutral pions in AuAu collisions for two centrality 
sets. The points are PHENIX data, histograms are 
the model calculations. 



Figure 4: The azimuthal two-particle correlation 
function for pp and for central AuAu collisions. 
The points are STAR data, dashed and solid his- 
tograms are the model calculations for pp and 
AuAu events respectively. 
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Figure 5: The nuclear modification factor, 
R AA (pr), for inclusive charged hadrons in central 
PbPb collisions triggered on jets with Ej l > 100 
GeV. The number of histogram entries and statisti- 
cal errors correspond to the estimated event rate for 
one month of LHC run. 
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Figure 6: The nuclear modification factor, 
R > AA (pr), for jets of cone size R = 0.5 in central 
PbPb collisions. The number of histogram entries 
and statistical errors correspond to the estimated 
event rate for one month of LHC run. 



4. Jet quenching at LHC 

The developed model was applied to analyze various novel features of jet quenching in heavy 
ion collisions at the LHC. Let us give a few examples of such jet observables. All calculations have 
been done for PbPb collisions at ^Ann = 5.5 TeV with PYQUEN energy loss model. The jet was 
defined on the generator level by a simple way, just collecting the energy around the direction of a 
leading particle inside a cone R = \/ Ar] 2 + A<p 2 = 0.5. The pseudorapidity cuts corresponding to 
the geometrical acceptance of CMS experiment were applied: | T) | < 3 for jets and neutral hadrons, 
1 7] | < 2.5 for charged hadrons and muons. 

4.1 Nuclear modification factors for jets 

The nuclear modification factor can be determined for jets by the same way as for inclusive 
hadron production. Since at the LHC no pp data will be available at ^/snn = 5.5 TeV at the time of 
the first PbPb data taking, particle spectra in pp collisions will be interpolated to this energy using 
perturbative QCD predictions constrained by the existing Tevatron data at 1.8 TeV and by the LHC 
results at 14 TeV. Another possibility to quantifies medium-modified particle spectra is to use the 
central to peripheral heavy ion collision ratio, Rcp{pt,ii), which does not require a pp reference, 
but has rather limited statistical reach of the peripheral data set. 

Figure || shows the /^-dependence of nuclear modification factor, R AA (pr), for inclusive 
charged hadrons in central PbPb events triggered on jets with Ej 1 > 100 GeV The number of 
entries and the statistical errors correspond to the estimated event rate for one month of LHC run 



and a nominal integrated luminosity of 0.5 nb [14]. The estimated suppression factor slightly 



increases with p T (> 20 GeV), from ~ 0.25 at p T ~ 20 GeV to ~ 0.35 at p T ~ 200 GeV. This 
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behaviour manifests the specific implementation of partonic energy loss in the model, rather weak 
energy dependence of loss and the shape of initial partem spectra. Without event triggering on 
high-is^ jet(s), the suppression is stronger (~ 0.15 at 20 GeV and slightly increasing with pj up to 
~ 0.3 at 200 GeV). 

A novel observable at the LHC will be the nuclear modification factor for hard jets, which 
can be reconstructed in high multiplicity environment with a good efficiency and low background 



starting from the energy ~ 50— 100 GeV [jl4j. Figure Q shows the pj -dependence of jet nuclear 



modification factor, R^ a (pt). The other conditions are the same as it was described above. The 
estimated suppression factor (due to partial gluon bremsstrahlung out of jet cone and collisional 
loss) is about 2 and almost independent on jet energy. The measured jet nuclear modification factor 
will be very sensitive to the fraction of partonic energy loss carried out of the jet cone. 

4.2 Medium-modified jet fragmentation function 

The "jet fragmentation function" (IFF), D(z), is defined as the probability for a given product 
of the jet fragmentation to carry a fraction z of the jet transverse energy. In nuclear (AA) interac- 
tions, the JFF for leading hadrons (i.e. the hadron carrying the largest fraction of the jet momentum) 



can be written as [ 15, 



d{p\?dydz< V ft) J (FT> y d(^fdy 



D ^)= j d{A) 2 dydz! ^^jj S[z-^ \ I I dtfffdy^tr, (4.1, 

Jet 
Z 'PT min 



where p\ = zp^ = z'pt is the transverse momentum of a leading hadron, z' is the hadron momen- 
turn fraction relative to the pj of the parent parton, pr T min is the minimum momentum threshold 
of observable jets, (dN AA ')/(d(p^ t ) 2 dy) and (dN AA )/(d(pj) 2 dydz') are the yields of^-typejets 
and hard hadrons, respectively. 

Figure [7] shows JFF's in central PbPb collisions with and without partonic energy loss for 
£^ et > 100 GeV. The number of entries and the statistical errors again correspond to the estimated 
event rate for one month of LHC run. Significant softening of the JFF (by a factor of ~ 4 and 
slightly increasing with z) is predicted. 

The medium-modified JFF is sensitive to a fraction e of partonic energy loss carried out of 
the jet cone, which is related also to the suppression of the absolute jet rates. Figure [8] shows 
the £-dependences of jet nuclear modification factor R S AA and ratio of JFF with energy loss to JFF 
without loss, D AA (z > zo) /£> pp (z > zo), for zo = 0.5 and 0.7 in central PbPb collisions [0 If e 
close to 0, then/? J ^ ~ 1 (there is no jet rate suppression), and JFF softening is maximal. Increasing 
e results in stronger jet rate suppression, but effect on JFF softening becomes smaller, especially 
for highest z (the ratio D AA /D PP can be even greater than 1 at large enough e and z values). The 
physical reason for the effect to be opposite in the jet suppression factor and the fragmentation 
function is it follows. Increasing e results in decreasing final jet transverse momentum (which is 



the denominator in definition of z in JFF (|4.1[)) without an influence on the numerator of z and, as a 
consequence, in reducing effect on JFF softening, while the integral jet suppression factor becomes 
larger. The crossing point between two effects is e ~ 0.3. 

Thus a novel concurrent study of the possible softening of the JFF and suppression of the ab- 
solute jet rates can be carried out in order to differentiate between various energy loss mechanisms. 
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Figure 7: Jet fragmentation function for leading 
hadrons (Ej > 100 GeV) in central PbPb collisions 
without (squares) and with (circles) partonic energy 
loss. The number of histogram entries and statisti- 
cal errors correspond to the estimated jet rate for 
one month of LHC run. 



Figure 8: Jet nuclear modification factor (solid 
curve) and ratio of JFF with energy loss to JFF 
without loss (z > 0.5 for dashed curve and z > 0.7 
for dash-dotted curve) as a function of the fraction 
e of partonic energy loss carried out of the jet cone. 



Strong JFF softening without substantial jet rate suppression would be an indication of small-angle 
gluon radiation dominating the medium-induced partonic energy loss. Increasing the contribution 
from wide-angle gluon radiation and collisional energy loss leads to jet rate suppression with less 
pronounced softening of the JFF. If, instead, the contribution of the "out-of-cone" jet energy loss is 
large enough, the jet rate suppression may be even more significant than the JFF softening. 

4.3 Jet azimuthal anisotropy 

The azimuthal anisotropy of particle spectrum is one of the most important tools to study 
properties of dense QCD-matter created in heavy ion collisions. It is usually characterized by the 
second coefficient of the Fourier expansion of particle azimuthal distribution, so called elliptic flow 
coefficient, vi. The momentum dependence of V2 for high-py hadrons, observed in semi-central 
AuAu collisions at RHIC, strongly supports the presence of rescattering and energy loss of hard 
partons in the azimuthally asymmetric volume of the nuclear reaction. A novel observable at the 
LHC will be the azimuthal anisotropy for hard jets (due to the part of partonic energy loss carried 
out of jet cone). 

The anisotropy of medium-induced partonic energy loss goes up with increasing collision 
impact parameter b, because the azimuthal asymmetry of the interaction volume gets stronger. 
However, the absolute value of the energy loss goes down with increasing b due to the reduced mean 
path length and the initial energy density. The non-uniform dependence of the loss on the parton 
azimuthal angle q> (with respect to the reaction plane) is then mapped onto the final parton spectra 
in semi-central collisions which are approximated well by the elliptic form [[T7|, [18|]. It results in the 
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elliptic anisotropy of observed high-p^ hadrons and hard jets. Figure |9] shows calculated impact 
parameter dependence of V2 coefficient for jets with Ej l > 100 GeV and for inclusive charged 
hadrons with pj > 20 GeV/c in PbPb events triggered on jets. The absolute values of vi for 
high-/?r hadrons is larger that one's for jets by a factor of ~ 2 — 3. However, the shape of b- 
dependence of v\ and v^ 1 is similar: it increases almost linearly with the growth of the impact 
parameter b and becomes a maximum at b ~ \ .6Ra (where Ra is the nucleus radius). After that, the 
vi coefficients drop rapidly with increasing b: this is the domain of impact parameter values, where 
the effect of decreasing energy loss due to a reducing effective transverse size of the dense zone 
and initial energy density of the medium is crucial and not compensated anymore by the stronger 
non-symmetry of the volume. 

4.4 Pj -imbalance in dimuon tagged jet events 

An important probe of medium-induced partonic energy loss in ultrarelativistic heavy ion col- 
lisions is production of a single jet opposite to a gauge boson such as a prompt y [|l^] or a y*/Z° 
decaying into dileptons [^, 21]. The advantage of such processes is that the mean (i.e. averaged 



over all events) initial transverse momentum of the hard jet equal to the mean initial/final transverse 
momentum of boson, and the energy lost by the parton in the QCD medium can be directly esti- 
mated from the observed pj -imbalance between the jet (or leading particle in a jet) and the lepton 
pair. 

In the 7+jet case the main problem arises from the jet pair production background when a 
leading 71° in the jet is misidentined as a photon. The "photon isolation" criteria usually used in pp 
collisions do not work with the same efficiency in high multiplicity heavy ion interactions. On the 
other hand, the production of jet tagged by dileptons is not affected significantly by backgrounds. 
The main background source - correlated semileptonic heavy quark decays - can be rejected using 



tracker information on the dilepton vertex position [14]. The moderate statistics, ~ 500 — 1000 



Z°/y*(— > jU + /i )+jet events per 1 month of LHC run with lead beams, are expected for the CMS 



geometrical acceptance and reasonable kinematic cuts Q14] , y.% . 

Figure pi] shows the difference between the transverse momentum of a jU + /i~ pair, p^, ^ , and 
five times the transverse energy of the leading particle in a jet (since the average fraction of the par- 
ent parton energy carried by a leading hadron at these energies is z ~ 0.2) for minimum bias PbPb 



collisions [22]. The process was simulated with CompHEP/PYTHIA generator package without 
and with partonic energy loss as obtained in the PYQUEN. The cuts p^>5 GeV/c, Pj ^ > 50 
GeV/c and Ef > 50 GeV, were applied. Despite the fact that the initial distribution is smeared 
and asymmetric due to initial-state gluon radiation, hadronization effects, etc., one can clearly see 
the additional smearing and the displaced mean and maximum values of the /^-imbalance due to 
partonic energy loss. The /^-imbalance between the pl + pi~ pair and a leading particle in a jet is 
directly related to the absolute value of partonic energy loss, and (unlike the /^-imbalance between 
the jU + jU~ pair and jet itself) almost insensitive to the form of the angular spectrum of the emitted 
gluons or to the experimental jet energy resolution [E2]. 
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Figure 9: The impact parameter dependence of el- 
liptic flow coefficients v^ 6 ' for jets with E S j l > 100 
GeV (black circles) and v\ for inclusive charged 
hadrons with pj > 20 GeV/c (open circles) in 
PbPb events triggered on jets. 
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Figure 10: The distribution of the difference be- 
tween the transverse momentum of a Z°/y* — > 
pair, pj ^ , and five times the transverse 
energy of the leading particle in a jet, 5£ , ^ ead , in 
minimum bias PbPb collisions with (dashed his- 
togram) and without (solid histogram) partonic en- 
ergy loss. 



5. Conclusions 



The method to simulate jet quenching in heavy ion collisions has been developed. The model 
is the fast Monte-Carlo tool implemented to modify a standard PYTHIA jet event. The full heavy 
ion event is obtained as a superposition of a soft hydro-type state and hard multi-jets. The model 
is capable of reproducing main features of the jet quenching pattern at RHIC: the pj dependence 
of the nuclear modification factor and the suppression of azimuthal back-to-back correlations. The 
model has been applied to analyze new features of jet quenching pattern at LHC energy: jet nuclear 
modification factor, jet fragmentation function, jet azimuthal anisotropy and dilepton-jet correla- 
tions. The further development of the model focusing on a more detailed description of low-pr 
particle production is in the progress. 
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